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Magnetoelectric (ME) materials are of fundamental interest and show broad potential for tech¬ 
nological applications. Commonly the dominant contribution to the ME response is the lattice- 
mediated one, which is proportional to both the Born electric charge Z® and its analogue, the 
dynamical magnetic charge Z™. Our previous study has shown that exchange striction acting on 
noncollinear spins induces much larger magnetic charges than those that depend on spin-orbit cou¬ 
pling. The hexagonal manganites i?Mn 03 and ferrites RFeOs {R = Sc, Y, In, Ho-Lu) exhibit strong 
couplings between electric, magnetic and structural degrees of freedom, with the transition-metal 
ions in the basal plane antiferromagnetically coupled through super-exchange so as to form a 120° 
noncollinear spin arrangement. Here we present a theoretical study of the magnetic charges, and of 
the spin-lattice and spin-electronic ME constants, in these hexagonal manganites and ferrites, clari¬ 
fying the conditions under which exchange striction leads to an enhanced Z™ values and anomalously 
large in-plane spin-lattice ME effects. 

PACS numbers: 75.85.+t,75.30.Et, 75.70.Tj, 75.47.Lx 


I. INTRODUCTION 

The cross-coupling between magnetic, electric, and 
elastic properties can lead to a plethora of novel and 
profound physical phenomena, with potentially broad 
and innovative applications. Magnetoelectric (ME) ef¬ 
fects are those in which the electric polarization P re¬ 
sponds to an applied magnetic field H, or magnetization 
M responds to an applied electric held S. The ME cou¬ 
pling (MEC) between magnetic and electric properties 
has motivated intense experimental and theoretical in¬ 
vestigations in bulk single crystals, thin hlms, composite 
layers, and organic-inorganic hybrid materials in recent 
years 

At the linear-response level, the linear MEC tensor a 
is dehned as 
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where indices /3 and v denote the Cartesian directions 
and fj,Q is the vacuum permeability. From a theoretical 
point of view, the linear ME effect can be decomposed 
into electronic (frozen-ion), ionic (lattice-mediated), and 
strain-mediated responses.^ Each term can be further 
subdivided into spin and orbital contributions based on 
the origin of the induced magnetization. As the orbital 
moment is usually strongly quenched on the transition- 
metal sites, most phenomenological and hrst-principles 
studies have focused on the spin-electronicji^ and the 
spin- lattic^™l contributions. The lattice response can 
be written, following Ref. [13 as 

= , ( 2 ) 

(sum over repeated indices implied), i.e., as a matrix 
product of the dynamical Born electric charge Z®, the 
inverse force-const ant matrix K~^, and the dynamical 


magnetic charge Z“, where m and n are composite la¬ 
bels for an atom and its displacement direction, fig is the 
unit cell volume. Note that Z“ is the magnetic analog 
of the dynamical Born charge, and is defined as 
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where Um is an internal displacement, Fm is an atomic 
force, and 77 is a homogeneous strain. In principle, Z™ 
has both spin and orbital parts, corresponding respec¬ 
tively to spin and orbital contributions to M,,, or Zee- 
man and p • A terms induced by iJi/, but we shall focus 
on the spin part in the following. Our previous first- 
principles study has shown that exchange striction acting 
on noncollinear spin structures induces much larger mag¬ 
netic charges than when Z“ is driven only by spin-orbit 
coupling (SOC). Therefore, exchange striction provides a 
promising mechanism for obtaining large MECs.^ 

The hexagonal manganites RMnOg and ferrites RFeOg 
(i? = Sc, Y, In, and Ho-Lu) form an interesting class 
of materials exhibiting strong couplings between electric, 
magnetic, and s tructu ral degrees of freedom.^^ A series 
of first-principle^i^lJ^l^ and phenomenologicaPi^ studies 
have greatly enhanced our understanding of the coupled 
properties. The ferroelectricity is induced by the struc¬ 
tural trimerization, and the direction of the spontaneous 
polarization is related to the trimerization pattern.^ An 
interesting “cloverleaf” pattern formed from interlock¬ 
ing domain walls between structural and ferroelectric do¬ 
mains has been found in hexagonal RMnOg and is now 
understood in terms of Landau theory.l^il^ Hexagonal 
RMnOg and RFeOg have rich magnetic phase diagrams 
and show considerable potential for manipulation and 
practical applicationsxhe magnetic order has two 
different origins, with the transition-metal Mn^+ or Fe^'*' 
sublattices ordering first, often followed by ordering of 
the rare-earth ions at lower temperature. The mag- 
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netic anisotropy is easy-plane and easy-axis for 3 d and 4 / 
spins respectively; the 3 d moments are antiferromagneti- 
cally coupled through superexchange so as to form a 120° 
noncollinear arrangement in the x-y plane, while the 4 / 
rare-earth moments are collinear along the hexagonal z 
axis. 


The low-temperature magnetic phases of i?Mn03 and 
i?Fe03 allow a linear ME effect to be present. The 
recently developed ME force microscopy technique has 
been used successfully to observe the ME domains in 
ErMn03.I^In that work, a large ME component Uzz ~ 13 
ps/m has been measured at 4 K, which is below the Mn^+ 
ordering temperature of 81 K but above the Er^+ or¬ 
dering temperature of 2 K. However, first-principles cal¬ 
culations predict that the SOC-induced spin-lattice azz 
arising from the Mn^+ order is 0 . 7 - 1.0 ps/m.^l^ This dis¬ 
crepancy suggests that the dominant ME effect in the 
hexagonal z direction is mediated by the Er^+ 4 / elec¬ 
trons in ErMn03. The in-plane ME effect, which has 
not been measured or calculated, has a different ori¬ 
gin. It is dominated by an exchange-striction mecha¬ 
nism on the Mn^+ sublattice, because the noncollinear 
spin pattern is sensitive to the lattice distortion. Thus, 
hexagonal i?Mn03 and i?Fe03 are good candidates to 
show exchange-striction enhanced magnetic charges and 
anomalously large spin-lattice MFCs. 

In this work, we use first-principles density-functional 
methods to study the magnetic charges and the spin- 
induced MECs arising from the 3 (i electrons in hexagonal 
HoMn03, ErMnOs, YbMn03, LuMn03, and LuFeOs. 
For the transverse magnetic charge components and 
MECs, we also provide a comparison between results in¬ 
duced solely by exchange striction and ones including 
SOC. Our results confirm that the exchange striction 
greatly enhances the in-plane magnetic charges, while the 
SOC contribution is minor for most components except 
on Mn atoms. However, the effect of SOC on the MECs 
is surprisingly large in many cases. This occurs because 
the exchange-striction contribution tends to be reduced 
by cancellations between modes, while the SOC contribu¬ 
tion is mainly associated with a few low-frequency modes. 
The in-plane ME responses are discussed case by case, 
and the conditions under which exchange striction leads 
to anomalously large in-plane spin-lattice MECs are clar¬ 
ified. 


The paper is organized as follows. In Sec. |H A| and 
|HB| we introduce the geometric structure and magnetic 
order of hexagonal i?Mn03 and i?Fe03. In Sec. m 
we analyze the tensor symmetries of the Born charges, 
magnetic charges and MECs in two different magnetic 
phases of i?Mn03 and i?Fe03. The computational de¬ 
tails are described in Sec. nmi The results and discussion 
of Born charges, magnetic charges and MECs in i?Mn03 
and LuFe03 are presented in Sec.|HI| We summarize and 


give our conclusions in Sec IV 



(b) 




FIG. 1. Structure of ferroelectric hexagonal JJMnOs (6 f.u. 
per primitive cell), (a) Side view from [110]. (b) Plan view 
from [001]; dashed (solid) triangle indicates three Mn®"*" con¬ 
nected via Opi to form a triangular sublattice at z = 0 
(2 = 1 / 2 ). 


II. PRELIMINARIES 
A. Hexagonal RMnOa 


Above the structural transition temperature Tc ^ 900 - 
1500 K, the hexagonal manganites i?Mn03 {R = Sc, Y, 
In, and Ho-Lu) are paraelectric insulators. The space 
group is POs/mmc with two formula units (f.u.) per 
primitive cell. Below Tc, the size mismatch between 
the small-radius ion and the large MnOs bipyra¬ 
mid leads to an inward tilting of the three corner-shared 
MnOs polyhedra and an associated “one-up/two-down” 
buckling of the R^'^ ion layer, as shown in Fig. The 
transition triples (“trimerizes”) the unit cell, lowers the 
structural symmetry to P63cm, and induces ferroelec- 
tricity. As the induced polarization is nonlinearly cou¬ 
pled to the trimeriz ation, these systems are improper 
ferroelectrics 

The Mn^+ magnetic order develops below the Neel 
temperature Tn of ~ 70 - 130 K. The in-plane Mn-O-Mn 
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FIG. 2. Magnetic phases of hexagonal RMnOs and i?Fe 03 . 
Mn^"*" ions form triangular sublattices at 2 = 0 (dash line) and 
2 = 1/2 (solid line), (a) A 2 phase with magnetic symmetry 
PGsc'm'; spins on a given Mn®"*" layer point all in or all out. 
(b) Ai phase with the magnetic symmetry PGscm, with Mn®^ 
spins pointing tangentially to form a vortex pattern. The Ai 
and A 2 phases differ by a 90° global rotation of the spins. 
The Bi and B 2 phases can be obtained from A 2 and Ai by 
reversing the spins on the dashed triangles. 


superexchange determines the noncollinear 120° antifer¬ 
romagnetic (AFM) order on the Mn^+ triangular lattice. 
On the other hand, the inter-plane Mn- 0 -i?- 0 -Mn ex¬ 
change, which is two orders of magnitude weaker than 
the in-plane exchange, modulates the relat ive sp in direc¬ 
tions between two consecutive Mn planes.l^^l^ At tem¬ 
peratures lower than ~ 5 . 5 K, the rare-earth ions with 
partially filled 4 / shells develop collinear spin order along 
the hexagonal 2 direction. For the Mn^+ order, there are 
four distinct magnetic phases, namely Ai (Pbacm), A2 
(Phsc'm'), Bi (POgCm'), and B2 (PGgc'm). The linear 
ME effect exists only in Ai and A2 phases. The Ai and 
A2 phases are shown in Fig. the Bi and B2 phases 
can be obtained from A2 and Ai by reversing the spins 
on the dashed triangles. From previous experiments, 
it is known that at zero temperature without a mag¬ 
netic field, HoMnOg is in the Ai phase, while ErMnOg, 
YbMnOg, and LuMnOg are not in either A phase. Under 
a weak magnetic field along the z direction, ErMnOg and 
YbMnOg undergo a transition into the A2 phase.l^^^^ 


B. Hexagonal AFeOg 

Epitaxially grown thin-film hexagonal i?FeOg has a 
similar structure as hexagonal i?MnOg, with improper 
ferroelectricity below ~ 1000 K. Replacing Mn^+ with 
Fe^"*" introduces larger spin moments and stronger super¬ 
exchange interactions in the basal plane. In a recent 
experiment, AFM order has been found to develop at 
Tn = 440 K followed by a spin-reorientation transition 
below Tr = 130 K in LuFeOg.^ It has also been con¬ 
firmed that below 5 K, the magnetic structure of LuFeOg 
is that of the A2 phase.^ 


C. Symmetry 

Our purpose is to understand the mechanisms that 
generate large magnetic charges that may in turn in- 
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FIG. 3. Symmetry patterns of Born charges, magnetic charges 
and ME tensors in AMnOg and AFeOg. (a) Tensor form of 
the ME conpling in the Ag phase, Born charges on Ai and 
Opi sites in either A phase, and magnetic charges on the 
same sites in the Ag phase, (b) Tensor form of the Born and 
magnetic charges on Ag and Op 2 sites in either A phase, (c) 
Tensor form of the Born charges on Mn, Fe, Oxi, and Ot 2 
sites lying on an My mirror plane in either A phase, and of 
the magnetic charges on the same sites in the Ag phase, (d) 
Tensor form of the ME coupling in the Ai phase, and of the 
magnetic charges on Ri and Opi sites in the Ai phase, (e) 
Tensor form of the magnetic charges on Mn, Fe, Oxi, and 
0x2 sites lying on an My mirror plane in the Ai phase. 


duce anomalously large spin-lattice MECs. Therefore, 
we focus on the Ai and Ag magnetic phases, shown in 
Fig. which allow a linear MEG to exist. ErMnOg, 
YbMnOg, and LuMnOg actually adopt other phases as 
their ground-state magnetic order at low temperature. 
Nevertheless, we include them for purposes of comparison 
when calculating the properties of the hexagonal AMnOg 
materials in the Ag phase. We also study LuFeOg in the 
Ag phase, and for HoMnOg we study both the Ag and 
Ag phases. 

The Ai and Ag phases have the same P6gcm structural 
symmetry, so the forms of the atomic Born charge tensors 
in the two phases are the same. The Born charges for 
Ri and Opi take the tensor form shown in Fig. |^a), 
while those of Ag and Opg have the symmetry pattern 
shown in Fig. [^b). For the Mn, Fe, Oti, and Ot2 sites 
lying on a vertical My mirror plane, the Born charges 
are as given in Fig. |^c); for the partner sites related by 
rotational symmetry, the tensors also need to be rotated 
accordingly. 

The symmetry forms of the atomic magnetic charge 
tensors can be derived from the on-site magnetic point 
symmetries. For the Ag phase, the magnetic space group 
is P6gcm and the magnetic charges of Ag and Opg take 
the forms given in Fig. [^d); those for Ag and Opg have 
the tensor symmetry shown in Fig. [^b); and for Mn, 
Fe, Oti, and Ot2 they can be written in the form of 
Fig. I^e). For the Ag phase, the magnetic group is 
P6gC^'; all the improper operators are associated with 
the time-reversal operation, so the magnetic charges have 
the same tensor forms as the Born charges. 

A symmetry analysis of the structure and the mag- 
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TABLE I. Atomic Born charge tensors Z® (in units of |e|) for 
LuMnOa and LuFeOa in the A 2 phase. TM = Mn or Fe. 



LuMnOs 

LuFeOa 


LuMnOa 

LuFeOa 

■^L(Lui) 

3.61 

3.79 

ZL(Oti) 

0.19 

0.11 

^L(Lui) 

4.12 

3.94 

ZUOti) 

-3.19 

-3.21 

•^SxjLua) 

3.66 

3.84 

ZUOt2) 

-1.90 

-2.15 

Zyx{ljU2) 

0.13 

0.15 

ZUOt2) 

-0.20 

-0.19 

^L(Lu2) 

3.96 

3.88 

Z^yiOT2) 

-1.85 

-2.13 

^S.(TM) 

3.17 

2.96 

ZL{Ot2) 

-0.18 

-0.11 

ZUTM) 

0.44 

0.21 

ZUOt2) 

-3.33 

-3.30 

Z^y (TM) 

3.26 

3.01 

ZUOpi) 

-3.00 

-2.40 

^L(TM) 

0.07 

-0.02 

zuopi) 

-1.54 

-1.61 

ZL(TM) 

3.95 

4.16 

zUOp2) 

-3.05 

-2.45 

ZUOti) 

-1.92 

-2.19 

^^.(0P2) 

-0.03 

-0.02 

ZUOti) 

0.25 

0.25 

^L(Op2) 

-1.43 

-1.52 

Z^y (Oti) 

-2.00 

-2.28 





netic space group identifies the phonon modes that cou¬ 
ple to the electromagnetic field. The infrared (IR)-active 
phonon modes that couple to the electric field are the 
longitudinal Ai modes and the transverse Ei modes, 

Tir = lOAi + 15Ei , (4) 

including the three acoustic modes. The magnetization 
is generated by phonon modes that couple to the mag¬ 
netic field. In the Ai phase, the magneto-active phonon 
modes are the longitudinal A 2 modes and the transverse 
El modes, 

= 5 A 2 + 15Ei , (5) 

where one pair of acoustic Ei modes are included. In the 
A 2 phase, on the other hand, the IR- and magneto-active 
phonon modes are identical, since the magnetic and Born 
charge tensors have the same form in this case. 

For the MFCs in the Ai phase, as the longitudinal 
IR-active and magneto-active modes are mutually exclu¬ 
sive, the ME tensor takes the form of Fig. [^d), which 
does not have a longitudinal ME component. For the A 2 
magnetic phase, the Ai and Ei modes are both IR-active 
and magneto-active, so that the ME tensor has both lon¬ 
gitudinal and transverse components and adopts the form 
shown in Fig. [^a). 

D. First-principles methods 

Our calculations are performed with plane-wave den¬ 
sity functional theory (DFT) implemented in VASI^ us¬ 
ing the generalized-gradient approximation parametrized 
by the Perdew-Burke-Ernzerhof functional .^3 The ionic 
core environment is simulated projector augmented 
wave (PAW) pseudopotentials,!^ and the 4/ electrons 
are placed in the PAW core. We use a Hubbard U = 
4.5 eV and J = 0.95 eV on the d orbitals of the Mn/Fe 
atoms, and the moment on the rare-earth ions are not 


TABLE 11. Eigenvalues of the force-constants matrix (eV/A^) 
for IR-active modes in LuMnOa and LuFeOa in the A 2 phase, 
excluding translational modes. 


Ai 

LuMnOa 

modes 

LuFeOa 

El modes 

LuMnOa LuFeOa 

4.24 

3.48 

3.32 

3.56 

7.44 

6.70 

4.68 

4.62 

8.74 

8.41 

6.73 

6.97 

11.51 

11.47 

7.35 

8.09 

14.01 

12.03 

8.63 

8.83 

15.60 

15.59 

9.56 

9.24 

22.66 

20.53 

11.36 

11.37 

25.87 

22.83 

12.46 

12.46 

35.82 

28.46 

13.02 

13.85 



14.09 

14.92 



16.49 

16.87 



17.37 

17.35 



23.36 

21.19 



37.75 

28.75 


considered.li^ The structures are fully relaxed in the 
DFT-|-tP^ calculations with their non-collinear spin ar¬ 
rangements in two cases, when SOC is present and when 
it is absent. In our noncollinear magnetization calcula¬ 
tion, a high cutoff energy 700 eV and a tight energy error 
threshold 1.0 x 10“® eV are necessary to get fully con¬ 
verged magnetic properties. The Born effective charge 
tensors and the F-point force-constant matrices are ob¬ 
tained using linear-response methods in the absence of 
SOC. The dynamical magnetic charges are computed by 
applying a uniform Zeeman fielcP^l to the crystal and 
computing the resulting forces. Polarization is calculated 
using the Berry phase formalism.!^ A 4x4x2 F-centered 
k-point mesh is used in the calculations. 


III. RESULTS AND DISCUSSION 
A. Born charge and force-constant matrix 

The / electrons are not included in our calculations for 
the hexagonal RMnOa class of materials, so the major 
differences between compounds result from the variation 
of the rare-earth radius; the trimerization tends to in¬ 
crease as the radius of the rare-earth element decreases. 
Because of the similarity in the geometric structures, the 
dielectric and phonon properties are almost identical in 
the RMnOa compounds, regardless of the magnetic or¬ 
dering. In Tables [T| and |TT| we list the Born charge tensors 
and the eigenvalues of the force-constant matrix for the 
IR-active modes of LuMnOa and LuFeOa. Only small 
differences are observed between LuMnOa and LuFeOa, 
reflecting the different transition-metal atom. The re¬ 
sults for the other i?MnOa compounds are quite similar 
to those of LuMnOa and are given for completeness in 
the Supplement. 
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TABLE III. Longitudinal magnetic charge components TABLE IV. Transverse magnetic charge components 

(10“®^b/A) of JJMnOs and LuFeOa in the A 2 phase. All (10 “^^b/A) of HoMnOa in the Ai phase, as computed in¬ 
components vanish in the absence of SOC. culding or excluding SOC. 



HoMnOs 

ErMnOs 

YbMnOs 

LuMnOg 

LuFeOg 


Total 

No SOC 


Total 

No SOC 

zZiRi) 

-50 

-53 

-53 

-67 

7 

^“(Hoi) 

-25 

-28 

^z”(Oti) 

-188 

-230 

Z7z{R2) 

14 

35 

24 

16 

7 

2”(Ho2) 

-15 

-18 

^((!r:(OT2) 

-57 

-67 

Z“(TM) 

-92 

-86 

-61 

-67 

9 

^“(Ho2) 

-1 

3 

^^1;(Ot2) 

-20 

-26 

ZTziTM) 

24 

1 

6 

25 

2 

^“.(Mn) 

92 

54 

^r,(OT2) 

-192 

-231 

ZS^z(Oti) 

-49 

-44 

-41 

-19 

23 

Z^iMn) 

-10 

2 

^”.(Opi) 

-483 

-551 

ZTz{Oti) 

99 

81 

53 

33 

22 

ZTyiMn) 

41 

48 

^.”(Op2) 

395 

461 

Z^z{Ot2) 

-7 

-12 

-12 

-12 

0 

Z’^xiOri) 

23 

28 

^”.(Op2) 

184 

253 

ZTz{Ot2) 

-119 

-94 

-64 

-49 

-25 

^“(Oti) 

-7 

-7 




^S(Opi) 

-276 

-257 

-230 

-190 

54 







ZTziOp2) 

141 

140 

125 

100 

-35 








B. Magnetization and magnetic charge 

In the A 2 phase, the trimerization induces not only 
an electric polarization, but also a weak ferromagnetism 
in the z direction arising from out-of-plane tilting of the 
Mn3+ spin moments induced by SOC. The net magneti¬ 
zations in the 30-atom unit cell for A 2 -phase HoMnOs, 
ErMnOg, YbMnOg, and LuMnOg are 0.309, 0.303, 0.292, 
and 0.268 /tb, respectively. These magnetic moments are 
found to depend almost linearly on the tilting angle of 
the MnOg bipyramids, which takes values of 5.03°, 5.07°, 
5.16°, and 5.21° respectively in these four compounds, 
but in any case the variation is not very large. In con¬ 
trast, the result for LuFeOg is -0.077 /xb, which is much 
smaller and of opposite sign compared with the i?MnOg 
materials. 

The magnetic charges defined in Eq. ([^ are more sen¬ 
sitive to the local environment, and now the differences 
between i?MnOg compounds are more significant. We 
divide the magnetic charge components into two groups 
that we label as “longitudinal” and “transverse” depend¬ 
ing on whether the coupling is to magnetic fields along 
the z direction or in the x-y plane respectively.' 

The longitudinal magnetic charge components are cal¬ 
culated with a magnetic field directed along z, which 
is roughly perpendicular to the spin directions. These 
components are only non-zero when SOC is considered. 
The scenario here is similar to the case of a transverse 
magnetic field {H^ or Hy) applied to Cr20g, since the 
magnetization is along the z axis for Cr20g. It is there¬ 
fore not surprising to find that the longitudinal magnetic 
charges of i?MnOg and LuFeOg in Table in are compa¬ 
rable to the SOC-induced transverse magnetic charges in 
Cr 20 g.^The longitudinal magnetic charges for Opi and 
Op 2 in LuFeOg are opposite to, and about three times 
smaller than, the ones in i?MnOg. These results explain 
the differences between i?MnOg and LuFeOg regarding 
the magnitude and the direction of the weak ferromag¬ 
netism, which is generated by trimerization distortions 
involving vertical displacements of Opi and Op 2 . 

For the response to transverse magnetic fields, both 


the field and the spins lie in the basal plane, so the dy¬ 
namical magnetic charges are driven by both SOC and 
exchange striction. As the exchange-striction strength 
can exceed that of the SOC by orders of magnitude in 
some materials, it is worthwhile to understand the rel¬ 
ative size of these two effects in i?MnOg and LuFeOg. 
In Tables [TV] and 0 we present the transverse magnetic 
charges induced with and without SOC in the Ai and 
A 2 phases. It is obvious that the SOC contributions 
are an order of magnitude smaller for many transverse 
components. Similarly, the magnetic charges induced by 
exchange striction are about ten times larger than the 
SOC-driven longitudinal ones in Table [llTl However, the 
SOC is crucial for the Mn atoms and it even reverses the 
signs of their transverse magnetic charges. 


C. Magnetoelectric effect 


We calculate the spin-lattice MFC from Eq. ([^ us¬ 
ing our computed Born charges, force-const ant matrices, 
and magnetic charges. The spin-electronic contributions 
are calculated based on the dP/dH version of Eq. Q 
with the lattice degrees of freedom frozen. We further 
subdivide the ME tensor components into longitudinal 
and transverse ones based on the direction of H relative 
to the hexagonal axis as before, so that the longitudinal 
(transverse) spin-lattice MEC is calculated using the lon¬ 
gitudinal (transverse) magnetic charge components. The 
MEC tensor elements allowed by symmetry are the lon¬ 
gitudinal azz and transverse Uxx = otyy ones in the A 2 
phase, and only the transverse ayx = —otxy components 
in the Ai phase. 

In the first part of Table 0!1 the spin-contributed 
longitudinal MFCs are shown for i?MnOg and LuFeOg 
in the A 2 phase. The MEC from the spin channel is 
dominated by the spin-lattice contribution. Although 
the longitudinal magnetic charges of LuFeOg are smaller 
than for i?MnOg, the spin-lattice MFCs \azz\ in i?MnOg 
and LuFeOg are similar, ^0.25ps/m. The results are 
comparable to those reported for the transverse MEC 
in Cr 20 j^ and for azz in ErMnOj^^ in previous first- 


principles calculations. In the second part of Table VI we 
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TABLE V. Transverse magnetic charge components Z™ (10 ^/tb/A) of ^RMnOa and LuFeOs in the A 2 phase, as computed 
including or excluding SOC. 



HoMnOs 

Total No SOC 

ErMnOs 

Total No SOC 

YbMnOs 

Total No SOC 

LuMnOs 

Total No SOC 

LuFeOs 

Total No SOC 

Z”(Ri) 

-23 

-24 

-21 

-22 

-37 

-40 

-42 

-35 

-36 

-52 

ZTx{R2) 

6 

-1 

6 

3 

12 

9 

14 

6 

15 

24 

Z“(R2) 

16 

18 

11 

12 

10 

10 

8 

7 

-9 

-11 

Z“(TM) 

-2 

10 

-7 

-10 

-16 

-21 

-11 

1 

-52 

-43 

Z“(TM) 

-42 

-24 

-38 

-22 

-25 

-34 

-31 

-17 

-102 

-95 

Z”(TM) 

-5 

46 

-7 

32 

-22 

27 

-32 

15 

-16 

-11 

Z”(Oti) 

5 

5 

6 

6 

12 

16 

14 

11 

0 

0 

Zz“(Oti) 

191 

221 

150 

154 

162 

178 

150 

122 

128 

105 

^.“(Oti) 

24 

23 

22 

22 

31 

33 

34 

25 

15 

11 

Z“(Ot2) 

20 

23 

16 

19 

19 

22 

17 

12 

25 

20 

Zz”(Ot2) 

195 

217 

140 

161 

173 

189 

166 

134 

130 

110 

(Ot2) 

-59 

-61 

-48 

-46 

-57 

-60 

-57 

-45 

-41 

-42 

Z“(Opi) 

-445 

-510 

-392 

-422 

-532 

-602 

-564 

-499 

-665 

-609 

Z“(Op2) 

241 

234 

215 

202 

298 

299 

316 

247 

388 

356 

Z”(Op2) 

-378 

-422 

-335 

-355 

-466 

-506 

-498 

-427 

-673 

-621 


show the spin-related transverse MECs for i?Mn 03 
and LuFeOa in the A 2 phase. The same information is 
presented in graphical form in Fig. 

It is clear from the comparison between the first and 
second parts of Table |VI| that the transverse spin-lattice 
MECs are one order of magnitude larger than the longi¬ 
tudinal ones, as a result of the exchange-striction mech¬ 
anism. Surprisingly, Fig. |^a) shows that the effect of 
SOC on the exchange striction is profound, even revers¬ 
ing the sign of the spin-lattice MECs in iJMnOa. This 
unusual behavior can be traced mainly to two observa¬ 
tions about the spin-lattice contributions from different 
IR-active modes in the RMnOa materials. Firstly, the 


TABLE VI. Computed MECs Qzz (longitudinal) and 
and ayx (transverse) for RMnOs and LuFeOs (ps/m). Spin- 
lattice, spin-electronic, and total spin couplings are given as 
computed with and without SOC. 


Spin-latt. Spin-elec. Total spin 

Total No SOC Total No SOC Total No SOC 




Oizz 

in 

As phase 



HoMnOs 

-0.27 

0 


0.06 

0 

-0.21 

0 

ErMnOs 

-0.26 

0 


0.05 

0 

-0.21 

0 

YbMnOs 

-0.25 

0 


0.06 

0 

-0.19 

0 

LuMnOs 

-0.19 

0 


0.00 

0 

-0.19 

0 

LuFeOs 

0.26 

0 


0.00 

0 

0.26 

0 




in 

As phase 



HoMnOs 

-0.99 

5.12 


4.10 

4.83 

3.11 

9.95 

ErMnOs 

-1.30 

2.40 


2.56 

3.72 

1.26 

6.12 

YbMnOs 

-2.52 

1.20 


3.72 

4.66 

1.20 

5.86 

LuMnOs 

-2.60 

1.31 


3.82 

3.50 

1.22 

4.81 

LuFeOs 

-2.20 

-1.57 


-0.79 - 

0.32 

-2.99 

-1.89 



0!.yx 

in 

Ai phase 



HoMnOs 

9.55 

4.88 


5.24 

5.35 

14.79 

10.23 


exchange-striction MEC is smaller than expected as a 
result of a large degree of cancellation between the con¬ 
tributions from different transverse IR-active modes. To 
illustrate this, the mode-by-mode contributions are pre¬ 
sented for a few selected cases in Table [Vlll Secondly, the 
softest modes are dominated by Mn displacements, pre¬ 
cisely those for which SOC has the largest effect on the 
Z™ values, even flipping the sign of some components. In 




HoMnOj ErMnOjYbMnOj LuMnOj LuFeO^ HoMnO^ 

(A,phase) 


FIG. 4. Transverse MECs for itMnOs and LuFeOs. a^x 
(ps/m) in the A 2 phase and ayx in the Ai phase, (a) Spin- 
lattice; (b) spin-electronic; and (c) total spin couplings. 
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TABLE VII. Transverse MEC contributions (ps/m) from IR- 
active modes for A 2 and Ai phases of HoMnOs and A 2 phase 
of LuMnOs. Results are given in ascending order of force- 
constant eigenvalues, which are reported in Table II of the 
Supplement. 


A 2 phase HoMnOa 
Total No SOC 

Ai phase HoMnOs 
Total No SOC 

A 2 phase LuFeOs 
Total No SOC 

0.01 

0.12 

0.25 

0.18 

0.28 

0.39 

-1.16 

2.62 

4.98 

2.36 

-0.54 

-0.50 

0.66 

2.32 

3.59 

2.37 

-1.31 

-1.22 

-0.51 

-0.35 

-0.32 

-0.48 

1.30 

1.23 

2.79 

3.13 

2.87 

3.33 

3.31 

3.12 

0.35 

0.21 

0.30 

0.30 

1.84 

1.73 

-1.88 

-1.85 

-1.35 

-1.90 

-4.43 

-4.11 

1.13 

1.25 

1.19 

1.38 

-2.59 

-2.25 

-2.96 

-3.07 

-2.70 

-3.40 

1.24 

1.13 

0.01 

0.13 

0.19 

0.06 

-1.48 

-1.27 

0.21 

0.24 

0.21 

0.26 

-0.15 

-0.14 

0.36 

0.40 

0.34 

0.42 

0.89 

0.83 

-0.03 

-0.03 

-0.03 

-0.04 

-0.62 

-0.55 

0.02 

0.01 

0.03 

0.03 

0.07 

0.03 


this way, it turns out that SOC can result in large rela¬ 
tive changes in the MEC results. In the case of LuFeOa, 
the SOC effect on the Z™ values is weak, even for Fe 
atoms. Thus, the MEC of LuFeOa does not change as 
dramatically as that of i?Mn 03 when SOC is included. 

From Fig. Qb) it can be seen that the spin-electronic 
contribution is not negligible in the transverse direction, 
and it counteracts the MEC from the spin-lattice channel 
in A 2 phase i?Mn 03 . The total transverse ME effect is 
summarized in Fig.j^c). Because of the large SOC effect 
and the cancellation between the lattice and electronic 
contributions, the total spin MEC axx is ~ 1.2 ps/m in 
most A 2 -phase i?Mn 03 compounds, except for HoMnOs. 
In HoMn 03 , the cancellation between the spin-lattice 
and the spin-electronic MFCs is the weakest of all the 
i?Mn 03 compounds, resulting in the largest total spin 
MEC of ^3.1 ps/m in the A 2 phase. In LuFe 03 , the 
spin-lattice and spin-electronic terms are all smaller than 
in i?Mn 03 . However, the cancellation induced by the 
SOC perturbation and the spin-electronic contribution is 
avoided, so that LuFe 03 has a large total spin MEC of 
~ —3 ps/m. 

We present the MFCs for HoMn 03 in the Ai phase in 


the last line of Table VI and in Fig. In principle the 
MFCs of HoMnOs in the Ai and A 2 phases should be 
the same without SOC, as the two phases only differ by a 
global spin rotation. This is approximately confirmed by 
a comparison of the corresponding entries for HoMn 03 in 
Table|Vll The ME contribution from exchange striction 
(i.e., without SOC) is ~ 5 ps/m for both the A 2 and Ai 
phases. However, when the effect of SOC is included, the 
spin-lattice contribution is strongly enhanced by another 
~ 5 ps/m. Furthermore, the spin-electronic MEC has the 
same sign as the spin-lattice one, which adds ~ 5 ps/m 
to the MEC. Therefore, the total spin MEC Uyx reaches 
~ 15 ps/m, and is the largest in all of the AMnOs and 
LuMn 03 materials we studied. 


IV. SUMMARY 

In summary, we have studied the spin-related magnetic 
charges and MFCs for HoMn 03 , ErMn 03 , YbMn 03 , 
LuMn 03 , and LuFe 03 using first-principles calculations. 
We confirm that the exchange striction acting on non- 
collinear spins induces much larger magnetic charges 
than does SOC acting alone. Nevertheless, the effect of 
SOC on the MFCs is surprisingly large, rivaling that of 
exchange striction in many cases. This occurs because 
the exchange-striction contribution tends to be reduced 
by cancellations between different IR-active modes, while 
the SOC contribution is mainly associated with just a few 
low-frequency modes with large Mn displacements. We 
also find that the AMnOs materials have spin-electronic 
MFCs comparable to the spin-lattice ones. Among the 
RMnOs and LuFe 03 materials we studied, we find that 
the Ai phase of HoMn 03 is the most promising ME ma¬ 
terial, with the largest MEC of ~ 15 ps/m. Extrapolating 
our conclusions to other hexagonal i?Mn 03 and i?Fe 03 
compounds that are not included in our calculations, we 
predict that the A 2 phase is more promising for the fer¬ 
rites, while the Ai phase has a stronger MEC for the 
manganites. 
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